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Abstract 

We develop a possible cosmology for a Universe in which there are n addi- 
tional spatial dimensions of variable scale, and an associated scalar field, the 
radion, which is distinct from the field responsible for inflation, the inflaton. 
Based on a particular ansatz for the effective potential for the inflaton and 
radion (which may emerge in string theory) , we show that the early expansion 
of the Universe may proceed in three stages. During the earliest phase, the ra- 
dion field becomes trapped at a value much smaller than the size of the extra 
dimensions today. Following this phase, the Universe expands exponentially, 
but with a Planck mass smaller than its present value. Because the Planck 
mass during inflation is small, we find that density fluctuations in agreement 
with observations can arise naturally. When inflation ends, the Universe re- 
heats, and the radion becomes free to expand once more. During the third 
phase the Universe is "radiation-dominated" and tends toward a fixed-point 
evolutionary model in which the radius of the extra dimension grows, but the 
temperature remains unchanged. Ultimately, the radius of the extra dimen- 
sions becomes trapped once again at its present value, and a short period of 
exponential expansion, which we identify with the electroweak phase transi- 
tion, ensues. Once this epoch is over, the Universe reheats to a temperature 
<m£if, the electroweak scale, and the mature Universe evolves according to 
standard cosmological models. We show that the present day energy density 
in radions can be smaller than the closure density of the Universe if the second 
inflationary epoch lasts ~ 8 e-foldings or more; the present-day radion mass 
turns out to be small (m ra dion^ eV, depending on parameters). We argue 
that although our model envisages considerable time evolution in the Planck 
mass, substantial spatial fluctuations in Newton's constant are not produced. 
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I. INTRODUCTION 



Recently it was suggested that the fundamental scale of gravity may be as low as TeV [|TJ. 
According to this idea, the observed weakness of gravity is associated with n new, relatively 
large spatial dimensions (compactified to a size ~ r ) in which only gravity can propagate. 
In this picture, all the standard model particles live in a set of branes with three extended 
space dimensions ("brane modes"), while gravitons live in the higher dimensional bulk of 
spacetime ("bulk modes"). This scenario turns out to be quite natural in (Type I) string 
theory. In this "brane world" picture |2j||, the standard model particles are open strings 
whose ends must end on the branes {e.g., stretched between branes), while gravitons are 
closed string states that can move away from the branes and into the bulk. The relation 
between today's Planck scale Mp\ = 1.2 x 10 19 GeV and the fundamental string scale m s is 
approximately given by 

Ml, ~ r<+ 2 r£ (1.1) 

Phenomenological and astrophysical constraints imply that m s may be as low as a few TeV, 
with n > 2 jiyij. In string/M theory, n < 7, and in the brane world, n = 2 is a reasonable 
choice ||. In any case, r must be fine-tuned to a very large value m s r ~ (M p i/m s ) 2 ^ n ^> 1; 
this fine-tuning problem is known as the radion problem. In this paper, we shall simply 
assume that the radius at is a stable minimum. However, as we shall see, this is not the 
end of the radion problem. One must still find a way for the radius to get to its final value 
without violating cosmological bounds. Some of the cosmological issues in the brane world 
have been discussed already pHTOf. 

The main concern of the present paper is the cosmology in this framework, during the 



epoch before big bang nucleosynthesis, especially inflation Obviously, in the brane 

world, the standard cosmological picture is altered dramatically. In this paper, we present a 
plausible cosmological scenario where a number of issues in the brane world, such as inflation, 
density perturbation, reheating, baryogenesis, as well as the radion problem, are addressed. 
Our scenario incorporates the brane inflation feature M and some of its extensions M, as 
well as some features of the rapid asymmetric inflation ||. The main goal here is to show 
that cosmology in the brane world is viable, and highlight some of the issues that we believe 
to be important. 

The reader may view this scenario as a search for viable potentials for the inflaton and 
the radion. The particular form we use has an effective potential for the radion field r and 
inflaton fields if) = (ipi, . . . , ipjsf) such that 

V(i/>, r) = Vb(V)[l + fi(r)] + / (r) + V^) (1.2) 

in the Jordan frame. (Appendix A gives some stringy justification for such a potential.) 
Here Vq(i/>) ~ m\ while Vi(ij)) ~ m^ w <C Vo, where mEw is the electroweak scale. Also 
the function /i(r) tends to force the radion to some value r\ when Vo ('»/') is large, whereas 
fo{r) is unimportant until r —>■ r , where r is the value of r today. Instead of choosing r$ 
to be the only minimum of fo{r), we choose /o( r ) to have multiple minima, with r just one 
of many possible minima. Our scenario utilizes an inflaton potential Vq(i^) due to brane 
separation (see Appendix B), and incorporates a radion potential fo(r) (and /i(r) as well) 
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with multiple minima. To simplify the problem, we might assume that Vo( , 0) depends only 
on one component of if) and Vi(ip) on the others (perhaps only one other); or, if) might only 
have one component with Vq and V\ both depending on this one component. The specific 



choices we examine for these potentials and functions are given in Eq. ( 2.27 ) for Vo(V0> ^q. 
( PH ) for fi(r), and Eqs. ( |2.68| ), ( |2.70| ) and ( |2.80| ) for fo(r). Here is a brief chronological 
description of the various phases of the scenario. 

Phase 0: The pre-infiationary phase: The key feature of this phase is that the radion 
field r is driven to a value 77 at which it becomes fixed, thus allowing the subsequent stage of 
standard inflation to take place. The piece of the potential (|1.2|) that achieves this trapping 
is the term Vo(ip)fi(r); rj will be a local minimum of fi(r). The initial conditions for this 
phase that we assume are that the radii of the extra dimensions begin at values "around" 
l/m s , where m s is the string scale (say, around 10 TeV). [By "around" we mean that values 
~ 1 — 100 times larger than m" 1 are not out of the question.] Such initial conditions are 
natural since in string theory, the only scale is the string scale, and thus all parameters 
should typically scale like m s unless there are good reasons (such as dynamical evolution) 



for other values. In § [IIB| we explore conditions under which the radion potential /i(r) 
achieves the fixing of the radion to some value r\. A specific choice of functional form 
of fi(r), which may or may not correspond to reality, is discussed further in Appendix [D]. 
Generally speaking, we think it likely that if fi(r) has numerous potential minima separated 
by some scale ~ m -1 , then the radion will become trapped at a minimum at fairly large 
mr (i.e. mr> a few), and that, if the radion settles to a potential minimum, it does so 
right away, without moving away from the potential well it starts in. The reason is that 
it is the Einstein- frame radion potential oc r~ 2n fi(r) that is relevant to the dynamics, not 
the Jordan- frame potential oc fi(r) [see Eqs. ( |C15| ) and fl2.10| ) below, and Appendix ^ for 
a discussion of the Jordan and Einstein frames]. It is plausible that fi(r) might have fixed 
amplitude of variation, so that the Einstein-frame potential will decrease oc r~ 2n at large r. 
Even under these assumptions about fi(r), the dual conditions that 77 should be "around" 
m" 1 and that mri> (a few), can be satisfied for values of m "around" m s . If fi(r) increases 
in amplitude at large r rapidly enough to overcome the factor r~ 2n , then it is possible that 
r actually increases somewhat from its original value during this pre-inflation era before 
settling into a minimum at 77. 

Phase I: Inflation at small Planck mass: Once the radion is fixed at the value 77, 
slow-roll inflation can take place. The value of the Planck mass during inflation Mpu is 
much smaller than today's Planck mass Mpi, as Mp n = (rj/r ) n Mp l , where r is the value 
of r today. In the brane world, brane inflation [0] is quite natural. (A brief review is given 
in Appendix B.) In the brane inflation scenario, when branes are separated by a distance d, 
an effective potential V(d) is generated by gravitational and other closed string exchanges 
between the branes. The distance d plays the role of the inflaton. On the other hand, d is 
related to the vacuum expectation value of a brane mode if) = m 2 d, so V(if)) is a function of a 
brane mode. In a particularly intriguing scenario, the electroweak Higgs field in the standard 
model plays the role of the inflaton ip. In the n = 2 case, the inflaton potential (schematically 
shown in Fig. 1), may be taken to have the following (oversimplified) qualitative form 

7o(V) + W) ~ - e~^ mi ) + Vi(V0 (1.3) 
with its minimum at if) ~ tuew ~100 GeV. 
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FIGURES 




For large d (but still much smaller than 77), V ~ m\ is very flat, so the inflaton slowly 
rolls down the potential towards small d. In standard cosmology, the number of e-foldings 
required to solve the flatness and the horizon problems is around 60, but for brane cosmology, 
the required number of e-foldings is different [see Eqs. ( |3.10| ) and ( |5.3|) 1, partly because the 
energy scale of inflation is m s <C vtigvt ~ 10 12 TeV, but also because the Planck mass Mpy 
during inflation is considerably smaller than today. [Fortuitously, the number of e-foldings 
required to solve the horizon and flatness problems may turn out to be about 60, see Eq. 

The amplitude of primordial density perturbations generated by quantum fluctuations in 
the inflaton field during this inflationary epoch is ~ mf/m/Afpy for the inflaton potential, 
Eq. (|1.3|). [See Eq. ( |2.38| ) in § [11 B[ ] Since Mpy ~ m s {m s ri) n l 2 , this amplitude is approx- 
imately ~ m s /mj(m s ri) n / 2 , and to achieve the measured amplitude ~ 10~ 5 of primordial 
perturbations |12|], we must require (m s rj) n ^ 2 ~ 10 5 m s /mj. For mj ~ m s , this would mean 
that m s ri ^> 1 (e.g. ~ 10 5 for n = 2), but larger mj is expected naively (see Appendix 0), 
and it is conceivable that m s ri ~ 1 — 100 (i.e. larger than one but not by a factor as large 
as 10 5 ). 

Phase II: Radion growth and radiation domination: At the end of inflation, we expect 
the brane to be heated (since the inflaton ip is a brane mode) while the bulk remains relatively 
cold. We expect the reheat temperature T r to be below the Hagedorn temperature, which 
is typically lower than m s (say by a factor of 3 to 10). Since m s » rriEWi this temperature 
T r should be above the electroweak phase transition critical temperature T c l . So the finite 
temperature effective potential V(ip,T r ) has a minimum at the origin, as indicated by the 
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dashed line in Fig. 1. The inflaton rolls past the minimum of the T = potential towards 
■0 = and is trapped there. The rate of cooling to the bulk (~ T" +7 /m™ +2 ) is very 
small, since T/m s is small. In this radiation dominated phase, the radion potential is 
negligible (which is not hard to arrange), while the inflaton field remains frozen. Under 
these conditions, we find that the cosmological model tends toward a fixed point in which 
the temperature remains nearly constant, while the radius r grows as a power of time [see 
Eqs. Q2.55D and ( p2 . 5 6| ) ] . While this powerlaw solution holds, the radion potential energy is 
unimportant compared with its kinetic energy, which, however, decreases with time (and 
hence with increasing radius). 

After sufficient time elapses, the kinetic energy of the radion field drops to a value 
comparable with the amplitude of the radion potential fo{r), and the growth of the radion, 
which is substantial up to this point, is halted. If fo(r) had only a single minimum, it would 
be a fantastic coincidence if (a) that minimum were precisely at the value r and (b) the 
powerlaw growth of r halted exactly when that minimum was encountered. Since powerlaw 
growth of r during the radiation dominated phase that follows inflation is generic in our 
picture, it seems that we must require that the radion effective potential fo(r) have multiple 
minima. Since it is inevitable that the radion kinetic energy becomes smaller than the height 
of its effective potential after some elapsed time, the radion must become trapped at one of 
the minima of fo(r) eventually. As an illustration, we consider a periodic radion potential; 
r is just one of the infinity of minima of this potential. That the Universe settles to r 
is a cosmological accident in our scenario, although it is natural for the radion to settle to 
some radius much larger than its value during inflation (and much larger than the string 
scale m' 1 ). Thus, the radion problem is not so severe in our picture, which accommodates 
growth of the radius of the extra dimensions to a large, but stable value very simply. What 
we do not explain is why r (and hence Newton's constant) has a particular value among 
the infinity of possibilities. However, we identify what conditions must be satisfied by the 
underlying physical theory for the Universe to settle at r [see Eqs. ( |2.75| ) and ( |2.76| )] . 

Phase III: Second inflationary era and electroweak phase transition: When T drops 
below T c l , the stable minimum of the inflaton will yield the spontaneous symmetry breaking 
of the electroweak model. It is reasonable to suppose that the electroweak phase transition 
is first-order. (This should be easy to arrange in models with multi-Higgs fields, e.g., the 
minimal supersymmetric standard model.) In this case, some super-cooling is expected and 
the actual phase transition happens during a period some time after T has dropped below 
T Cj i. In the mean time, we expect considerable dilution of the radion energy density as well 
as the bulk energy density. Indeed, we show that requiring the radion density at present 
not to exceed the critical density for a flat Universe constrains the number of e-foldings of 
this inflationary era [see Eqs. ( |3.12|) and ( |5.2|) 1. [Associated with the radion density today 
would be small-amplitude oscillations of Newton's constant at a high frequency ||13|| , |f[4}| ; see 
Eqs. ( p.88| ) and §|V[] This can be achieved by a short period of inflation during the super- 
cooling period, followed by either prompt or delayed reheating. The actual electroweak 
phase transition then takes place with the presence of nucleation bubbles. This allows the 
electroweak phase transition to complete and baryogenesis during this period can happen 
more or less as in the standard scenario [15j . Alternatively, baryogenesis can happen via the 
Affleck-Dine mechanism or some other mechanism |16|| . In fact, it may take place before the 
end of the second inflationary era. 
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The final reheat temperature can be around a few GeV, maybe even close to the elec- 
troweak scale m EW , if m s is large enough, and still avoid excessive cooling to the bulk, 
which would overproduce Kaluza-Klein (KK) modes, whose energy density could overdose 
the Universe and ruin the success of Big Bang nucleosynthesis. It is also high enough to 
provide the "initial conditions" of the hot Big Bang before Big Bang nucleosynthesis. 

The scenario is summarized in Fig. 2. 

The basic plan of this paper is the following. In §pj we present our cosmological scenario; 
II A| gives some useful background (some of which is also found in Appendix |Cj), § |II B| treats 



the pre-inflationary phase, §11 (J treats inflation at small Planck mass, and § |11D| treats the 
phase during which the radius of the extra dimensions grows from r\ to r . Some constraints 
on our model are gathered in § |TTT| . Density fluctuations during the epoch when the radion 
grows are discussed briefly in §[TV| The results are discussed in §|V|. Some additional details 
about our cosmological model are contained in various Appendices. 



II. MATCHING PHASES 



A. Setup 

The starting point for our analysis is the following low energy action, which is valid 
when the lengthscales over which all fields vary are much larger than the size of the extra 
dimensions: 

R 



S = j d A X\J —g 

+ Srest[e '^g a p, Xrest 



(2.1) 



This action is derived from the higher dimensional description in Appendix |C]. The action 
is written in the Einstein frame, g a p is the Einstein frame metric, and 



9a/3 = e-^gofi (2.2) 

is the physical, Jordan frame metric. The quantity G is the usual 3-dimensional Newton's 
constant, and \i is a mass of order the Planck mass M P i = 
by 



ch/G = 1.22 x 10 19 GeV given 



M, 



pi 



In + 2 
327m 



(2.3) 



where n is the number of extra dimensions. The field <3> is the canonically normalized radion 
field, related to the radius r of the extra dimensions by 



r exp 



(2.4) 



where r is the equilibrium radius of the extra dimensions today. The fields if) = (tpi, . . . , ip^) 
are inflaton fields. The quantity V(ip, $) is the Jordan-frame potential for the radion 
and inflaton discussed in the Introduction and in Appendix [A] below. Finally the action 
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Srest[g a f3, Xrest] is the action of the remaining matter fields Xrest, which in our analysis below 
we will treat as a fluid. 

As discussed above, we assume that the Jordan-frame effective potential V(ip, $) = 
V(iJ),r) is of the form [cf. Eq. Q] 



V{ij>,r) = Vo(if>)[l + fi(r)] + / (r) + V x {$). (2.5) 

Here we might assume that Vo depends only on one component of ij) and V\ on the other 
components (perhaps only one other). The function fi(r) tends to force the radion to 
some value ri while Vo is large, whereas fo(r) is unimportant until r — > ro, its value today. 
Following Ref. @], we assume that Vb(?/>) ~ as -0 ~ and that Vb(i/>) asymptotes to 
a constant value, Vo, exponentially with some mass scale(s) mi. The potential Vi(i/') is 
assumed to have a minimum at nonzero if), and a value Vi = Vi(0) at if) — which satisfies 
Vi « V^o- 

We will do all our calculations in the Einstein frame. Assuming zero spatial curvature, 
the metric of the cosmological background in the Jordan frame can be written in the form 

d s 2 = -dt 2 + a 2 (t)dx-dx, (2.6) 

and the corresponding metric in the Einstein frame is 

ds 2 = e*^ \-dt 2 + a 2 (t)dx-dx\ = -dt 2 + a 2 (t)dx-dx, (2.7) 

where 

di = exp($/2/i) dt d(i) = exp($/2/i) a(t). (2.8) 



It is important to keep the relations ( |2.8| ) in mind, since proper time is not the same in the 
Einstein and Jordan frames, nor is the scale factor. In terms of r, the scalings are 

dt = {r/r ) n / 2 dt a(t) = (r/r ) n/2 a(t). (2.9) 

Thus, as r changes, the relative rates of advance of time and scale factor differ in the two 
frames. 

We treat the last term in the action (|2.1|) as a fluid with Jordan-frame density p and 
pressure p. Then, the cosmological equations of motion that follow from the action ( [2 . 1\ ) 
follow from the general equations of motion given in Appendix [C], and are given by [] 



H 2 ^ 



3M& 



2 2 



(2.10) 



These evolution equations do not include any coupling between the inflaton and the radiation, 
which would be necessary to describe reheating. If we add the standard type of phenomenological 
terms to achieve this in the Jordan frame, we obtain after transforming to the Einstein frame that 
one should add a term —F(t/j)i(j'e~^^ 2 ^ to the right hand side of the inflaton equation ( p,12| ), and 
a term F(i/j)^' 2 e^^ 2 ^ to the right hand side of Eq. ( |2.13D , where F(ip) can be any function. 
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d 

$" + 3H& + — 
<9$ 



-2$ 



^V(ip,$) 



2/i 



2/j, 



(2-11) 



3# — 



+ e-*/^V(V» ) *) = 0. 



and 



p +3(p + p) 



2// 



0. 



(2.12) 



(2.13) 



In these equations primes denote derivatives with respect to Einstein-frame proper time t 
[Eq. ( |2.6| ) above], and H is the Einstein-frame Hubble parameter H = a! /a. 

In the next few subsections, we shall substitute the potential (|1.2| ) into the evolution 
equations fl2.10|) - ( p.!3|) , and solve for approximate solutions in the four different phases of 
cosmological evolution discussed in the Introduction. 



B. Phase 0: Radion to Its First Equilibrium 



During Phase 0, the radion evolves to some size r\ where it remains pinned during 
inflation (phase I). This pinning happens by some time to when the scale factor is a . The 
end of phase signals the onset of the first inflationary phase, phase I 

To see how the pinning might come about, let us assume that we can neglect the terms 
fo(r) and Vi(if>), and that there is no energy density except what is due to if) and the radion. 
Furthermore, let us assume that the kinetic energy of the inflaton is negligible, and that 
Vo(V') ~ Vq. Then the Friedmann equations (|2.10| ) and (|2.11|) together with the potential 
(O) simplify to 



H 2 



57T 



3M|! 



V exp(-2$//i)[l + /i(r)] + -(# 



>\'2 



2V r exp(-2$//i) 



Mr) 



2 

r dfijr) 
2n dr 



We can render these equations non-dimensional by letting 

C = exp[2($ - $,)//!] 
and defining a new time variable r by 



dr 



3MfJ *M-*i/»)dt, 



(2.14) 



(2.15) 



(2.16) 



where $j is the initial value of at the beginning of Phase 0. In terms of these new 
variables, we find (/' = df/df in these equations) 



i + Mf) l /Q 2 
C 2vlc 



Mr) 



r dfj{r) 



2n dr 



(2-17) 



S 



Here y = a/a^ with fij the initial value of the Einstein frame scale factor, r = r i (, l l 2n ) with 
Ti the initial value of the radius of the extra dimensions, and 

4S " (2.18) 



n + 2 

It is easy to see that if fi(r) = 0, the solution of Eq. (|2.17|) tends to 

y oc f 1 / 3 C oc f^/ 3 , (2.19) 

according to which r grows without bound, and the radion kinetic energy dominates the 
energy density of the Universe, but there is no inflation [ 17j . The approach to this asymptotic 



solution may be very slow: for example, for n = 2 we have ( oc f 4 ^, so the vacuum energy 
density declines oc f 4 ^, which is only a bit faster than the rate of decline of the radion 
kinetic energy, 1/f 2 . Nevertheless, it is noteworthy that without the radion potential, r 
would grow to infinity in this phase. 

To halt the growth of r, the radion potential must be capable of trapping the radion 
field. From Eq. ( [2.17| ), we see that this is only possible if the condition 

can be satisfied. If we consider the choice fi(r) = (mr) K for example, where m is a mass 
parameter and k > 0, then it is clear that Eq. (|2.20|) will only have real solutions for k > 2n 
in which case the radion could settle to a value mr = (n/2n — l) -1 . (Analogously shifted 
minima have been discussed by e.g. Steinhardt & Will f!4|| .) However, such steeply growing 
functions f(r) could prevent r from growing to a large value later on. 
Instead, we shall consider the possibility that 

/i(r) = aF(r), (2.21) 

where a is a dimensionless amplitude factor, and F(r) has multiple minima, separated by 
a characteristic scale ~ m -1 , with "potential barriers" |-F(r)| ~ 1. A specific example (but 
not unique or required) is F(r) = 1 — cos mr, for which Eq. (|2.20|) becomes 



amr sin mr 

l+a(l- cos mr) = 0. (2.22) 

2n 

Since cos mr < 1 and sin mr < 1, there are no solutions to this equation unless amr/2n > 1, 
or mr > 2n/a. Thus, unless a is large (which we consider unlikely), the radion will only 
settle into minima at relatively large values of mr, if at all. This conclusion ought to hold 
for other choices of F{r) with similar qualitative properties. If, for example, m ~ m s , then 
we conclude that the radion will only settle on values larger than the string scale, which, in 
fact, is required for consistency of our entire picture. (Remember that, for example, r must 
exceed the brane separation.) 

For fi(r) of this general type, it seems likely that the radion must settle into its first 
minimum, the one nearest its value at the onset of Phase 0, if it settles to a minimum at 
all. The reason is that the height of the Einstein frame effective potential decreases with 
increasing r, so that if the radion acquires sufficient kinetic energy to roll over the first 
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barrier it encounters, it should be able to overcome all subsequent barriers. (Remember 



that in the asymptotic solutions to Eq. [|2.17|| the energy density of the Universe becomes 



dominated by radion kinetic energy as time progresses.) We explore the particular example 
F(r) = 1 — cos mr in some detail in Appendix [D]. It is also important to note, though, that 
it is possible to imagine choices for fi(r) that undo the decrease factor in amplitude cx r~ 2n 
in Eq. (|2.14j) at large r. For such potentials, it might be possible for r to evolve considerably 
before settling to a potential minimum. 

C. Phase I: Inflation at Small Planck Mass 

During Phase I, the Universe inflates at a fixed radion radius, r\. This fixes the Planck 
mass to be 

Ml, = {n/r ) n M* u (2.23) 



and the expansion rate is 



H 2_ (*\ 2 _ Sttt^Vq _ 87rr%V 



a 



in the Einstein frame, where prime denotes differentiation with respect to t. The Jordan 
frame expansion rate is simply 



Hl= -) = ^ L , (2-25) 



where dot denotes differentiation with respect to t; this is just as in the usual general 
relativity but with a different Planck mass. The scale factor in this phase is 

^ = exp[Hi(i - to)] = — = exp[Hj(t - to)]. (2.26) 
ao ao 

Since the radius of the extra dimensions is frozen in phase I, it is equally easy to use the 
Jordan or Einstein frame descriptions. (The same will not be true for subsequent phases.) 

Presuming that a single inflaton field is important in this phase, with an effective poten- 
tial 

V (if>) = V [l - exp(-^/m I )], (2.27) 
the inflaton equation of motion is 

$ + 3Hii> = -— exp(-VVmi). (2.28) 
mi 

The evolution of the inflaton proceeds in two subphases. During the first subphase, we have 
the usual slow-rolling approximation, 

4, « _ Wz^ll , (2 . 29) 
3Hjmi 
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whose solution is 



Mt -to) , , , , Mi 



PLI 



exp(^/mi) = exp(^oMi) Q 2 „ = exp(^ /mi) - hsH^t - t ). (2.30) 

Officii j oTTTflj 

This approximate solution holds as long as the two conditions 



Mg 1|I exp(-^/m I ) ^ ^ Mg M exp(-2^/m I ) i g 

247rm 2 V> 2^0 ~ 487rm 2 



are satisfied. Assuming that Mpy 3> mi, which emerges naturally later, the first of these 
two conditions fails first, when 

M 2 

expWmO ~ (2-32) 



If the initial value exp^o/mi) of exp^/mi) is far larger than the limit ( p. 32 ), then the time 



required for the inflaton field to reach this magnitude is extremely large, given byf] 



Hj{t BI - t ) ~ exp(VoM), (2.33) 

m pi,i 

where £ sr is the time at the end of slow-roll, which implies many e-foldings during inflation. 
When slow-rolling ends, a second sub-phase of inflaton evolution begins. Since the kinetic 
energy of the inflaton at the beginning of this subphase is only ~ (127rm 2 /Mp n )Vo <C Vo, the 
inflaton moves on approximately a "zero energy solution with negligible damping" . That is, 
it satisfies the equation 



i> « -v / 2V r exp(-^/2mi), (2.34) 

which has the solution 



exp(ip/2mi) = exp(ip sr /2mi) — 



\ 77^0 (* ~ *sr) = exp(^ sr /2mi) - J f PI>I - t s 
\ 2mf V Idtt mi 



(2.35) 



where exp(?/> sr /2mi) ~ Mp^i / mi\/247r . The time remaining for ip — > is not large: Hj(t cn d — 
t ST ) ~ 1, where t en d signifies the end of this inflationary epoch, and henceforth we do not 
distinguish between the two times t sr and £ en d. 

Using our approximate expression for Hi(t sr — t ) « Hj(t ew i — t Q ), we rewrite the slow 
rolling solution as 

M 2 

exp(^/mi) « -^4#,(t cnd - t). (2.36) 



2 The subscript "sr" stands for the end of slow rolling. 
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From this and the slow rolling approximation it follows that 



^ - -T^r, (2.37) 



The primordial density perturbation amplitude is then (18 



p ip miMpy 

where Nk is the number of e-foldings that remain between horizon crossing for a scale of 
comoving length ~ h~ x and the end of Phase I. We note that the spectrum of inhomogeneities 
implied by Eq. (|2.38| ) is insensitive to k, in agreement with observations [n|. If the effective 



potential as a function of interbrane separation d is proportional to 1 — exp(— m^d), and we 
set ip — m 2 d, then m^d = mdi/j/m 2 ., which implies m\ = m 2 /md. Then the amplitude of the 

density fluctuations is proportional to \J^md/m 2 Mpi i i, which, for Vo ~ m*, is ~ m^/Mpy. 
Note also that since the radion is trapped, fluctuations in r are suppressed [see Eq. ( |D9| ) of 
Appendix Q. 



D. Phase II: Radiation Domination 

At the end of inflation, t = i\ and a — a\\ the Universe reheats to a temperature eV^ A 
where e < 1 depends on how efficiently the kinetic energy of the inflaton is thermalized once 
ij> ->0. 

Reheating will alter the effective potential so that there is a minimum at ip = 0, provided 
that the critical temperature T Cj i for the phase transition (presumed first order) connected 
with the potential Vi(i/>) is small compared with the reheating temperature. In this case, 
there will be a nonzero vacuum energy Vx, but as long as the temperature remains above 
T c i, the Universe remains radiation-dominated, and if) is pinned at zero. 

We assume that the radion potential becomes negligible when this happens. Let us 
explore the growth in the radion field that ensues. 

It is most convenient to work in the Einstein frame. This phase of evolution will divide 
into three subphases. During the first subphase, when radiation dominates, the Friedmann 
equations are [cf. Eqs. fl2~T0l ) and fl213|) above with p = p/3] 

H 2 = J^pexp(-2$/7t) pa 4 exp(-2$//i) = constant, (2.39) 

where we recall that p is the Jordan frame energy density i.e. p ~ T 4 . Thus, the scale factor 
in Einstein frame is simply 

/ i \ 1/2 

a = &i(j-) > (2-40) 

just as it would be in constant-M P1 cosmology. However, to find out how fast the temperature 
decreases, we also need to know how the radion field evolves. 
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Under the assumption that the radion potential is negligible, we find from Eqs. ([L2] 
( p. 11| ) and Q2.27 ) the evolution equation 



2V 

$" + 3#$' = — exp(-2$/V). 



We make the change of variables 



Then Eq. ( |2.41| ) is equivalent to 



u = exp(2$/,u). 



u'\ . „ /3 u'\ , AVi 
u" + 3H )u' = u" + W = -4, 



(2.41) 



(2.42) 



2.43) 



u J \2t u / jj, 

where H = l/2t, appropriate for this subphase, has been used. During the first subphase 



3 u' 
2t u 



(2.44) 



In that case, under the assumption that v! = and u = u x at t = t x , we find 

5/r ' ' 

which can be integrated to yield 



(2.45) 



" = " 1 + ^( t - 5f " + l^ 

This solution holds as long as Eq. fl2.44j ) is true, a condition that fails when 

2tu' ^ 16V^t 2 /15/i 2 
3m ~ m + 4V r 1 £ 2 /5 / u 2 ~ ' 



^2.46) 



or 



(2.47) 



;2.48) 



(In getting this condition, we have assumed that the first subphase ends at a time ^> ti.) 
Since u = exp(2$//x) is still only a factor of two or so different from Ui, it follows that 



P 



piafexp(2$//i) piViiju p^ij 



a 4 exp(2$i//i) /x 2 -u 2 p?u\ 
where in the first equality we used Eqs. ( [2.40 ), ( 2.42 ) and ( 2.4S ). UsingF] 



^2.49) 



'Recall that M 2 n = M^(ri/r ) n = M 2 1V ^T. 
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2 3Mg LI Mg lv ^IT 
t l ~ ~ . 



Pi 



and t\ = y/ultf, we find 



p ~ Vi 



[2.50) 



(2.51) 



since /i ~ Mp\. Thus, the first subphase ends when the energy density in radiation becomes 
comparable to the vacuum energy in the inflaton. The value of the Planck mass only changes 
by a factor of order unity during this regime. 

During the second subphase, the radion potential is still unimportant, but the radion 
kinetic energy is not. Thus, the radion evolves according to Eq. (|2.43| ), but the expansion 
rate in the Einstein frame is given by 



H 



in 



3M^ 



P | ViW , P 2 
u u 8 



l\2 



2^1 



The evolution of the inflaton field is governed by the equation 



ijj" + 3Hij)' 



dip 



(2.52) 



(2.53) 



We assume that the Universe remains hot enough that the inflaton is trapped in a symmetric 
phase, at fixed fixed vacuum energy density V± during this subphase. At least at first, 
the kinetic energy of the inflaton field will be unimportant, and we can approximate the 
expansion rate by 



H 2 



S7T 



3M2> 



u 



8m 2 



(2.54) 



It is easy to see that at the start of this subphase, the three contributions to the energy 
density are comparable to one another. Moreover, there is a simple, powerlaw solution to 
the fully nonlinear problem defined by eqs. ( |2.43| ) and ( 2.54|) . For this solution, a oc 
and 



u 



u 
u 



2 
T 



(2.55) 



since pa A /u =constant, p =constant. From Eq. ( ^.54[ ) we find a consistency condition 

p 9n — 6 



n + 2 



(2.56) 



which gives p/V\ = 3 for n = 2, for example; thus, the temperature remains ~ Vl'\ and 
could be comfortably above the critical temperature for the inflaton in this regime.0 Within 
this solution, it also follows that, up to a possible additive constant, 



4 If we let a(t) = do(i)[l + cn(t)] and u(t) = uo(t)[l + T](t)], where ao(i) and uo(t) are the powerlaw 
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and therefore, as a function of t, the radion expands according to 



T//U4 

u = , (2.58) 

and grows without any expansion of the Universe at all in the Jordan frame. 

It is worth investigating the meaning of this powerlaw solution further. Clearly, a so- 
lution in which the temperature remains constant is not expanding in the Jordan frame at 
all. Such a solution ought to apply only in a limiting sense. That is, the correct solution 
might approach this one asymptotically, at late times. To see if this happens, we con- 
sider the numerical solution eqs. ( |2.43|) and ( |2.54| ), with the energy-conservation condition 



pa 4 /n=constant. We arbitrarily choose u = Ui and p = pi at some initial time tj, when 
a = cij. If we define 

y = ^ and C = -, (2.59) 

and define a dimensionless Einstein-frame time variable by 

df = (SnpjZM^UiY^di, (2.60) 
then we find the two coupled equations 

y) 2v\() C V 4 

C+(—-%)C' = w, (2.61) 



y C 



where 



solution of eqs. ( EU3D and (gjg) , it is easy to show that 



v"+ 3 -i + ^ 2 + 6 4- 

2t t 2 t 



, a ( 3n — 2 
a + 



r V 4n 



n + 2 
24n 



2t 



these coupled perturbation equations have powerlaw solutions oc t s where 

0. 



. t \ ( 2 s 3n — 2 

Thus, there are no growing perturbations, and the powerlaw, fixed point solution to eqs. ( p. 43 ) 
and ( 2.54 ) is stable. 
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48U and v = *. (2.62) 



n + 2 pi 
We can also evaluate the Jordan frame scale factor from 

a = w 1/4 a, (2.63) 

and can find the time elapsed in the Jordan frame using 

dt = u 1/4 dt, (2.64) 

which we rescale by defining tj = (8n pi/ SM^y/uiY^t, to get t'j = C -1 / 4 . Eqs. ( [2.61D are to 
be solved with initial conditions ( — y — 1 at f = ff, we can choose t« arbitrarily (although 
we expect it to be ~ 0.5 if we want a solution that has y — > at f — > 0). The solutions 
depend only on the two parameters, v and v ; in terms of these, the powerlaw solution found 
previously becomes 

Numerical evaluations for n = 2 show that, although the solutions oscillate slightly at late 
times, they approach the simple solution found above to high accuracy. From examining 
the output, it appears that the Jordan frame scale factor actually does not remain precisely 
constant at late times, but increases and even decreases slightly (by less than 10%) as time 
progresses. 

In the asymptotic regime of the second subphase, the radion grows according to 

AV f 2 

u = ( r /r ) 2n = -4-, (2.66) 
p 2 



so u — > 1, or r — > ro, at a time 

p 



*r — 



(2.67) 



— 1/2" 

if this solution continues to hold. Note that t r ^ ro ~ u l ti, where tj marks the onset of this 
subphase (or the end of the first subphase; e.g. Eq. |[2.48|| ); since U\ <S 1, t r ^ ro 3> ij. Q 

In order for the radion to become pinned at r = ro we need a coincidence to happen: near 
the time t r ^ ro , the radion potential itself must begin to play a central role in the evolution 



5 As this is a pretty odd solution, let us also consider an alternative, that the radion evolves like 
a free field after the second subphase begins, and its kinetic energy dominates the energy density 
of the Universe. In this case = ^af/a 3 , and 

which implies the solution (a/aj) 3 = t\/12vr < I>J/Mpi. Using this solution, we find that <J>' = 
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of the field. Only the radion's potential can make it settle into a minimum, rather than 
rolling forever to ever increasing radius. Indeed, what we want is for the effective potential 
of the radion to have many possible minima, so that the value it settles into eventually is 
determined by this coincidence. 

To understand the settling process better, we need to incorporate the radion potential 
term fo(r) of the Jordan-frame potential (|1.2p in out analyses. If Vbuik(^) is the radion 
potential in 4 + n dimensions, then, after integrating over the n extra dimensions, the 
corresponding Jordan frame potential is 

fo(r) = r n V hulk (r) (2.68) 

[see Appendix |C|, and the Einstein frame potential Ve obtained after conformal transfor- 
mation is 

V E = r 2 n r- n V hulk (r) = r w exp(-$///)y bulk (r) = (2.69) 

Ju 



Suppose that 

Vbuik(r) = U hulk F(r), (2.70) 

where [/bulk is a constant and F(r) is dimensionless. We assume that F(r) may undulate 
up and down, but with a characteristic amplitude |-F(r)| ~ 1; thus the scale of the radion 
potential is determined by £/buik- Since F(r) is dimensionless, it must contain mass scales; 
these are reflected in the magnitude(s) of the derivative(s) of the potential. We assume that 
F(r) may have multiple minima (an infinite number in the model considered below.) 

One condition for the radion to be able to settle into one of the minima of its potential is 
that the Einstein frame kinetic energy density fall below the Einstein frame potential energy 
amplitude. Since the kinetic energy is 

fi 2 fu'\ 2 [i 2 



8 \uj 2t 2 



and the potential energy amplitude is 

r-Q^buik _ r%U bu i k n 



u 



(2.72) 



2t\JV 1 

during the powerlaw regime, the two become comparable at at an Einstein frame time 



Mpi/i\/12n and consequently 

r=r \v 

or r oc i 1 /^ oc for n = 2. For this solution to hold true, all other contributions to the energy 
density must decline more rapidly than ($') 2 / 2 « i~ 2 . But it is easy to see that /9exp(— 2<J>//x) oc 
a -4 oc t -4 / 3 according to this solution, so it must not be valid. 
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t 



mm 



ro^buik 

This time ought to be smaller than or comparable to the time at which u 



require 



which implies 



ro^buik 



r^buik 



// 



2Vl 



(2.73) 
1 i.e. we must 

(2.74) 
(2.75) 



Since, presumably, i/buik an d Vi are determined by fundamental physics, this relationship 



may be taken to determine r^. Moreover, since we know that M\ 



pi 



m 



2+n„n 
f 0> 



we find 



2V 1 m 



2+n 



or 



bulk 



bulk 



2V l m 2+n 



(2.76) 



If we assume that V± ~ ^ew^ w h ere m £w is the electroweak unification scale, and m s > 
rriEWi then it is clear that U^^/mf^ 71 ~ m^^/M^mf <C 1. In getting these estimates, 
we have presumed that F(r) takes on a typical value for r ~ ro i.e. we are excluding the 
possibility that, for example \F(r ~ ro)\ <C 1, which would alter the above estimates. This 
amounts to assuming that whatever the mass scales appear in F(r) are generally of order 
Tq 1 or larger. 

To investigate the settling process in more detail, we need the equations of motion, which 
can be obtained from Eqs. ( j2.10 ) and ( 2.1 If ) . The resulting equations are 



H 2 



57T 



u 



3H 



u'\u' 



u u 
AV X 2TV l ^u 



U J u 



8u 2 

r dF(r) 
n dr 



F(r) 



r=r 



u l/2r, 



(2.77) 



These equations are augmented by the conservation condition pcfi/u =constant. Also we 
have defined 



T 



O^bulk 



(2.78) 



which we expect to be ~ 2. If we nondimensionalize as before we find 



>\ 2 



2v\( 



t\ 2 



C y 4 



vTy/v ~ 

C l/2 



F( ri ( 



l/2n\ 



vv 



1/2 



r dF(r) 
n dr 



F(r) 



[2.79) 



Notice that in this form of the equations, T only appears in the combination T 4 = Ty 7 ^, 
and since we expect Mj < 1, this parameter is small, implying that deviations from the 
powerlaw solution only appear at late times, as we have already concluded. 
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These equations contain three parameters explicitly: n (or u), v = V±/pi, and Tj. In ad- 
dition, they contain one (or more) parameters implicitly, because of the mass scales implicit 
in F(r). For example, if 

F(r) = l-cos(m r r), (2.80) 

so that there is only one mass scale, m r , there is an additional nondimensional parameter 
Hi = m r ri. For this form of F(r), the evolution equations are 

V'V 1 {C'\\v 1 »Ti 1/2 



f » + (^-|)c' = ^-^i C " 2 



/iiC 1/2n sin(/iiC 1/2n ^ 



- 1 + cos(/^( 



l/2n> 



f2.81} 



n 

In these units, substantial deviations from the scaling equations are expected after 



r - T,-y=±£ (2.82) 

at which time ( = 4/T 2 , provided that fii > (T i /2) 1 / n . 

Fig. 2 shows the results of numerically integrating the dimensionless evolution equations, 
eqs. (pU), for (v,T f ,^) = (10~ 6 , 10~ 8 , 10~ 3 ) (Fig. 2a) and (v, T is fa) = (10~ 4 , 10" 6 , 10" 2 ) 
(Fig. 2b), respectively, with n = 2 in both cases. The numerical results show clearly that 
after a long period of powerlaw expansion (in close agreement with the fixed point solution 
found above), ( levels off, although in both cases, the time at which this happens is a 
bit later than our back-of-envelope estimate, so that ( is systematically larger than 4/T 2 
asymptotically. For Tj = 10~ 6 , we would estimate ( = 4 x 10 12 asymptotically, whereas the 
numerical result is 9.75 x 10 13 , a factor of about 25 larger; for Tj = 10~ 8 we would estimate 
( = 4 x 10 16 , as opposed to the numerical 4 x 10 17 , about a factor of 10 larger. (The 
discrepancy in r is smaller, since ( oc r 2n = r 4 for n — 2.) These results can be explained 
if the time to asymptote is a factor of three to five larger than our simple estimate. By 
inspecting the figures, we can see that the time at which ( levels off is about five times larger 
than the analytic estimate, ~ 4.08 x 10 7 , for (v, T i; /ij) = (10~ 4 , 10~ 6 , 10~ 2 ), and about three 
times larger than the analytic estimate, ~ 4.08 x 10 10 , for (v, Tj,/ij) = (10~ 6 , 10" 8 , 10~ 3 ). 

Once the radion field begins to settle into a minimum around tq, the value of the Planck 
mass zeros in on its present value. Q Once this happens, the temperature of the Universe 



6 It turns out that for a quadratic potential, there is also an exact solution in the radiation 
dominated era if we ignore the source term. In this case, the equations are completely linear, and 
if we take a potential like ^m 2 (^ — ^o) 2 the solution is 

$ - $ 0- , (2-83) 

the most important feature of this equation is that the amplitude falls like (mi)~ 3 / 4 (i.e. oc a~ 3 / 2 ) 
at late times. 
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can begin to fall once more, and ultimately it must drop below T Cj i. When this happens, 
the inflaton fields once again are free to roll, and move toward their minimum at nonzero 
vacuum expectation values. The amplitude of any residual oscillations in the radion field 
will then redshift away exponentially, until the inflaton kinetic energy is thermalized in a 
second phase of reheating. It is therefore necessary that once the temperature of the Uni- 
verse becomes constant during the radiation-dominated era of radion growth, the constant 
temperature must be above T C) i. Moreover, we need to require that there is enough inflation 
after the temperature falls below T c l for the amplitude of oscillations in radius to drop to 
an acceptable level. 





2 4 6 8 10 

log 10 (H i t) 

FIG. 2. Numerical solutions for the cosmological model during Phase II for two different sets 
of parameters. The long-dashed line is Cw> the solid li ne is the dot-dashed line is a(t), and 
the dotted line is t{t). Both solutions are characterized by an initial adjustment period during 
which the radiation energy density, which is proportional to aj 4 , drops, followed by a protracted 
period of powerlaw growth of the radius of the extra dimensions (at fixed temperature), which 
terminates when the radion potential traps the field at a minimum. Exponential inflation begins 
once the radion becomes trapped. 



For F(r) = 1 — cos m r r, it is easy to see that the minima of Ve(3>) are at m r r = 1i\k r , 
where k r is an integer. At minima, 

<PV E {<$>) 47r 2 fcM*W r y 
and, at the minimum corresponding to r = r , 
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bulk 



H 2 n 2 



Thus, near r = r (or $ = 0) 



2n 2 k 2 r%U hn i k 2 
Ve($) w ~7~~ <3r 



fx 2 n 2 



and the mass of the radion is 

_47c 2 k 2 r%U huik 8n 2 k 2 ro Vi 25Gn :! ^V- 



To 



radion 



\i 2 n 2 



/j 2 n 2 n(n + 2)M } 



2 ' 
PI 



^2.85) 



^2.86) 



(2.87) 



numerically, we find (recall Eq. ||2.75|| ) 



m rad ion ~ 7.3 x 10 3 eV 



"r 



n(n + 2) VlTeV^ 



^2.88) 



Remarkably, the mass of the radion that emerges is much smaller than any other character- 
istic mass scale in the problem, unless k ro ^> 1 and/or V\ ^> 1 TeV 4 . For other choices of 
F(r) we would have 



d 2 V E ($) _ r%U hulk fr 



d$ 2 

with the definition 



H 2 n 2 \ r 



-n(n — l)F(r) + r 



,d 2 F(r) 



dr 2 



(2nk r 



—n(n — l)F(r) + r 



,d 2 F(r) 



dr 2 



(2.89) 



(2.90) 



this becomes identical to the formula for the special case F(r) = 1 — cos m r r. 



III. EXPANSION FACTORS AND THE RADION DENSITY 

Now that we have a complete account of the various phases of expansion in our proposed 
cosmological model, we can gather the results to calculate the factors by which the Universe 
has expanded between various interesting epochs and the present. We shall work in the 
Jordan frame, for the most part, and define the cosmological scale factor a (to) = 1 at the 
present day, t = t - The value of the Hubble constant today is H = 100/io kms" 1 Mpc -1 , 
the CMBR temperature today is T = 2.7° K, and the corresponding critical density is 
p Cj0 = 3H 2 M 2 JS7r « 8.01 x 10- 47 /i 2 GeV 4 . 

Before moving on to our more complicated cosmological model, it is useful to review the 
situation in convention cosmology with a fixed Planck mass and a single inflationary era. Let 
a exp be the value of the scale factor at the end of the period of exponential expansion, and 
a r h the value of the scale factor at the end of the reheating phase that follows exponential 
expansion. If T r h is the reheating temperature, then 

a rh ~ 1/3 , (3.1) 
<->rh 1 rh 
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where the dimensionless factors <5>o and $rh count particle states in thermodynamic equilib- 
rium at present and at a r h, respectively. [] If the energy density during inflation is pv, and 
the inflaton potential is harmonic near the minimum attained at the end of inflation, then 

-^o— ~^UJ' (3 ' 2) 



where £ r h is another dimensionless factor that counts the contributions of various particle 
states to the total energy density at temperature T r h- Q Defining 

2 c T 4 \ 1/4 



we find a exp /a rh ~ e% 3 , or 

4/3 9 l/3 T 1/3-1/3^1/4 , 2 T 4 x l/4 

°-'~^r = ^2^W) • (3 - 4) 

If we assume that the present day Hubble scale Hq 1 passed outside the horizon during 
inflation, then the scale factor of the Universe at that time was anubbie = (Pc,o/ Pv) 1 ^ 2 -, an d 
the ratio 



= exp 30.8 + 0.25 In ( T ^ TJ ] + ln' rh rh 



l/3 c l/3^1/4 2rpA 
dexp £ r h O t-rh 11 1 PV 



& Hubble S^[ 3 V 15/?c,0 



ITeVV V 



(3.5) 



For pv ~ (10 15 GeV) 4 , Eq. ( |3.5| ) implies about 60 e- foldings between auubbie and a e:C p, the 
familiar value, but generally the number of e-foldings depends on details of the inflationary 
model. 

In the cosmological model developed above, there are two periods of exponential inflation 
that occur at different values of the Planck mass. The second inflationary epoch, and the 
subsequent reheating, occur at the end of Phase II, at which time the Planck mass has 
settled to its present value. Thus, we can apply the same reasoning to this epoch as was 
developed in the preceding paragraph and we find that the scale factor at the end of the 
period of inflation that concludes Phase II is 

1/3 q l/3pl/4 4 1/4 



7 Recall that it is entropy that is conserved during adiabatic expansion. Note that S r h may be 
considerable e.g. >10. 

8 When only relativistic particles are present, and the equation of state depends only on temper- 
ature, £ r h = 2>S r h/A. 
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in the Jordan frame, where a r h,2i £rh,2 

and Slh 2 have the same meanings as the analogous 
symbols introduced for conventional inflation and reheating, but apply to the end of Phase 
II only. If we assume that this second inflation led to an increase in scale by a factor JF 2 , 
then the scale factor just before inflation began was 

1/3 <jl/3c.l/4 .^4. 1/4 
dex P ,2 e rh,2°0 C rh,2 ( nl 0\ /, 7 ^ 

this is also the value of the scale factor before the radion field began its powerlaw growth 
during Phase II. 

Proceeding backward in time still further, we encounter the first subphase of Phase 
II, during which the energy density of the Universe declined from its value just after the 
reheating at the end of Phase I, it 2 £ r h^T!f h 1 /15, to ~ V\, the corresponding increase in scale 
was a factor ~ (7r 2 ^ r / l l T r 4 h l /15Vi) 1//4 , so the scale factor at the end of the reheating that 
terminated Phase I was 

1/3 q l/3 F l/4 T 1/3 q l/3 p l/4 , 2 T 4 x l/4 

e rh,2^0 Wfe,2 J _ Crh^O ^r/1,2 / ^ ^ \ 7 / Q o\ 

a ^ ~ s i/3 f .1/4 T - 01/3 - — vT^v^ ) ' { } 

where e r h t i = (ix 2 £ rhy iT^ h l / IhVo) 1 / A . Consequently, the value of the cosmological scale factor 
at the end of the exponential expansion in Phase I is 

1/3 1/3 ol/S.l/4 2 T 4 X 1/4 

aez P ,i ~ e rM a rM —j^- I — rr I . (3.9) 

Although the derivation of a exPt i is more complicated than the derivation of a exPt 2 (or its 
conventional equivalent, a exp ), notice that it does not depend explicitly on the different 
Planck scales that arise in our inflationary model. 

The dependence on Planck scales enters when we reconsider the relationship between 
a exP) i and a Hubble- We have assumed that the present Hubble scale, and all other macroscopic 
scales relevant to the development of large scale structure, crossed the horizon during Phase 
I. As a result, a Hu bbie = (Pc,o/Vo) 1/2 (M Pli i/M P i), and 



1/3 1/3-1/3-1/4,, 2 „4 T -v 
dexp,l e rh,l € rh,2°0 C rh,2 M Pl( ^ 1 V 



a H ubbie SJ^JWpi,! V 15p£„ 



1/4 



exp 



1/3 1/3 c l/3 c l/4 



30.8 + 0.25 lnf + lnf -^-) + lnf € J±^lh^ ^ 



(3.10) 



For our cosmology, the number of e-foldings between horizon crossing and the end of expo- 
nential expansion during Phase I depends on numerous uncertain parameters, principally V 
and the combination Mpi/J^Mpy. The requirement that a exPt i > aHubbie is one constraint 
on our cosmological model. 

We can obtain a separate constraint by requiring that the energy density in radions 
today does not overfill the Universe. At the end of the subphase of powerlaw growth of the 
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radion during Phase II, the energy density in radions is ~ V\\ during the ensuing exponential 
expansion it drops to ~ V\jT\. The energy density in radions drops by an additional factor 

off Q>exp,2/ Q>rh,2)^ ~ ^Jh 2 by the time reheating is complete, to a value ~ ix 2 Erh^T^h 2/15^2 
a factor ~ T\ smaller than the energy density in relativistic matter at the end of reheating. 
Between T r h,2 and T , the radion density drops by a factor of SqTq / Srh^T^h 2; so that 



Pradfi ~ ^3 

J -0<-V/i,2~ / ~ 2 



is the density of radions today. Comparing with the closure density implies a radion density 
parameter 



thus, fi rod)0 ^ 1 as long as T rh)2 /J^ ^ IO^oS^/So^m) eV 



IV. FLUCTUATIONS 

Fluctuations about the smooth cosmological background alter the form of the Jordan 
frame line element from Eq. ( [2.6|) to 



ds 2 = -dt 2 + a 2 (t)(5ij + hi^dx^x 1 ; (4.1) 



the corresponding Einstein frame metric is ds 2 = <is 2 y w(x, t), where w(x, t) = 
exp[2<3>(x, t)/fi]. Assume that u(x, t) = tto(£)[l + r?(x,t)], where |?y(x, i)| <C 1. Then af- 
ter making an appropriate infinitesimal coordinate transformation, we get 

ds 2 = -di 2 + a 2 (t)(8ij + h ij )dx i dx i (4.2) 

i.e. the metric in the Einstein frame can be reduced to synchronous form. As always, the 
necessary infinitesimal coordinate transformation is not unique; there is still gauge freedom 



even when the metric is reduced to the form of Eq. (|4.2j ) | 19| , |20 |. 

The perturbed Ricci tensor corresponding to Eq. ( |4.2j ) can be found in Ref. [I9|. For 
compressional perturbations around the powerlaw background solution for Phase II, the 
perturbation equations can be boiled down to two gauge-independent equations, 



y, W H V 2 <p H fn + 2\ip 



2t 3a 2 V 24n J t 

(4.3) 





f9n - 6\ 


V- 






V 32n ; 




V 16n J 



V 



where is the (Einstein-frame) gauge-independent metric potential introduced in Ref. [20 
and tp — 7/ + 77/i, coupled to three gauge-dependent equations, 
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„ 37/ 77 V 2 77 ft' 

77 + — + ^ 

%' r /Qn _ fix A 

ft + 



2t £ 2 a 2 t 
9n-6^ A /n + 2\ 77 /n + 2W 



8n / £2 V 8n / 1 2 V 2 / £ 2 



(tf)' + = 0, (4-4) 



where h — ha. If po(t) and «o(*) denote the background solutions given by eqs. ( |2.56|) and 

^4^4 = ^[l+A(x,f)], (4.5) 
u(x,t) Uo{t) 

and, if the Einstein-frame three- velocity is U = W(x, t), 

x(1>i)= WM. (4.6) 



Eqs. ( [4.4| ) and the second of eqs. ( |4.3| ) can be combined to yield the gauge-independent 



equation 

„„ I3ip" /37t2-2\^' /9n-2\^ 3V 2 ^ V 2 ^' 6V 2 ^ 4V 2 0' H ,„ ^ 
^" + -3-+ — r— 7- + -j— = —s£ T^T->- ( 4 - 7 ) 



2t \ An J i \ An J t 3 to 2 a 2 t 2 a 2 to 2 

On scales larger than the Einstein-frame horizon scale, H" 1 , a complete solution may be 
obtained by coupling Eq. (|4.7|) to the first of eqs. (|4.3|). 

At the end of Phase I, there are no fluctuations in the radion field, so ip = on all scales 
of interest today (i.e. well outside the horizon). (This is because the effective mass of the 
radion field once in a minimum is generally larger than the cosmological expansion rate; 
for a particular example, see Appendix |B|.) There are fluctuations in 0# on these scales, 
and, from Eq. (|4.7|) , these tend to generate perturbations in the radion field. However, the 
driving terms are very small on large scales: for comoving wavenumber k, they are of order 
^(fin/Pa 2 <C (pn/t 4 - Consequently, if) ~ k 2 i(p H /a 2 <C 4>ii/t, and the large-scale fluctuations 
in the radion field generated during Phase II are far smaller than 0# (which sets the scale 
of the density fluctuations on such scales after they re-enter the horizon). Moreover, the 
source term in the equation for §n is negligible on these scales, since ip/t <C (pn/t 2 , and 4>h 
remains constant on scales larger than the horizon during Phase II. 



V. DISCUSSION 

In the preceding section, we developed a new picture for cosmology in the brane world. 
Our cosmological model is based on a specific form of the effective potential, Eq. fll.2| ), which, 
although admittedly somewhat complicated, allows the size of the compact dimensions of 
the Universe to evolve to its present value from a different, but much smaller, fixed value at 
early times. In the specific scenario we have unfolded, the resulting evolution of the Universe 
divides naturally into four different phases, the last of which can be called the "standard Big 
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Bang cosmology" that follows the electroweak phase transition and proceeds to the present 
day, with the radius of the extra dimensions fixed at its present value, r , and hence the 
Planck mass fixed at G~ 1/2 = 1.22 x 10 19 GeV. 

The other three phases represent the evolution of the radius of the extra dimensions to 
that value from a considerably smaller one. While we do not claim that the scenario we have 
developed for this evolution is unique, it does have some features that are attractive. The 
first phase, Phase of §|11 Bl is relatively brief; the radion settles into a potential minimum 



at T\ during this phase. We have shown that this process is not entirely guaranteed to take 
place, but may be rather likely in scenarios where the effective potential for the radion has 
multiple (or an infinite number of) minima. 



Phase sets the stage for Phase I of Q during which the Universe inflates at a fixed 



Planck mass M Pi) i = Mpi(ri/r ) n/2 < M P1 . We assume that this is the main inflationary 
phase undergone by the expanding Universe, so that macroscopic co moving scales on which 
large scale structure develops all passed outside the horizon during Phase 0. The density 
perturbation amplitude produced by quantum fluctuations in the inflaton field(s) if) during 
Phase is estimated in Eq. ( |2.38| ). For an inflaton effective potential proportional to m A s [l — 
exp(— rridd)], where d is the inter-brane separation (as discussed in |2T| and Appendix B), 



we estimate that the primordial density fluctuation amplitude is ~ (m^/Mpy)./^ for a mode 
with comoving wavenumber k, with the number of e-foldings between horizon crossing 
for that mode and the end of exponential inflation during Phase 0. (In Appendix B, we 
argue for = iurr, the mass of the RR mode.) Nominally, we would expect m d < m s , 
leading to a density perturbation amplitude ^m s /Mpi t j, which would be woefully small for 
m s ~ TeV if Mpy = Mpi. An attractive feature of our scenario is that it allows Mpy <C Mpi. 
Turning the argument around, observations of large scale temperature fluctuations in the 
cosmic microwave background radiation [^] require m^/Mpy ~ 10~ 5 , so the radius of the 
extra dimensions during Phase I satisfies the constraint 

riW 2 MpiT , n5 m d 10 _n/ m d \ 



rj M) 



pi 



Mpi VI TeV/ v ; 



The radius of the extra dimensions during Phase I was considerably smaller than today 
if ma ~ 1 TeV and n < 7. One of the principal motivating factors behind our cosmo- 
logical model is the realization that the amplitude of primordial density perturbations is 
proportional to Mp^, and that perturbations at an acceptable amplitude are only possible 
if M P y < Mpi. 

Expansion from to tq occurred during Phase II, which is mainly radiation-dominated 
following the reheating that terminated Phase I. In our model, once the branes come to 
overlap, Vo(ij>) — > 0, which frees the radion to expand once more since the product term 
Vq(iI>) j 'j(r) in the effective potential ( |1.2|) is no longer active. The inflaton potential is then 
dominated by Vi(i{)), and we assumed that the initial reheating was sufficient to trap the 
Universe at small ifr at first, at a minimum with nonzero V\(ij)). Phase II naturally divides 
into three sub-phases, which was discussed in detail in §|TTJ. During the first subphase, the 
Universe expands at r ~ rj, until it cools sufficiently that an approximate equilibrium is 
attained, with comparable energies in radiation, radion kinetic energy, and vacuum energy 
density. Once this happens, a new phase of powerlaw expansion of the radius of the extra 
dimensions ensues at virtually fixed radiation temperature; see eqs. (|2.55| ) and ( |2.56|) . This 
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subphase ends when the radion becomes trapped in one of the many (or infinite) minima 
of its effective potential, fo(r) (see Eq. [L2| ); we assume that this minimum is at r , and 



eqs. ( 2.75|) and ( |2.76| ) estimate the radion vacuum energy density in the bulk required for 



this to be true. The third subphase of Phase II is the phase transition associated with the 
inflaton potential Vx(ip) in Eq. ( |1.2| ). During this phase, the Universe expands exponentially 
by an additional factor JF 2 , and reheats, finally, to a temperature T r h$- Conventional, non- 
inflationary Big Bang cosmology commences at this point. 

Another important constraint on our cosmological model is the requirement that the 
present day radion energy density does not dominate the total energy density of the Universe. 
Because the radion is trapped in a potential minimum toward the end of Phase II, it behaves, 
much like the axion, as a massive, cold dark matter particle; the effective radion mass is 
estimated in Eq. ( |2.88| ), and may be ~ 1 eV typically.0 Just after the reheating that ends 
the short inflationary period during Phase II, the energy density in radions is smaller than 
the energy density of the products of reheating by a factor « JF^ 3 . Requiring that radions 
not dominate the mass density of the Universe today implies, by Eq. ( |3.12|) , 

— (i^) 1/3 (lf) 1/3 ' 

where T r h 2 is the temperature of the Universe after this last reheating episode, and the 



remaining factors are ~ 1 in general; see §111 for details. Fj Thus, if the second inflationary 



epoch comprised more than about eight e-foldings, the present day density in radions would 
be negligible, but our model cannot be consistent with fewer than eight e-foldings, which 
would result in an overdense Universe dominated by radions. We note that Newton's con- 
stant of gravitation actually oscillates in this model at a frequency m radion , but with a very 
small amplitude, SG /G ~ ^l a \ H /m Iadion (e.g. jll). 

Combining eqs. Q3.10D , (|5.1| ) and ( |5.2| ), we can constrain the number of e-foldings that 



take place during Phase I between the time when the present-day Hubble length passed 
outside the horizon and the end of exponential expansion: 



— ~ exp 

(^Hubble 



/ritevx / luuGev y/* | / 

V m d ) V T rh , 2 J \ 



1/3 1/3 n l/3 



48.7 + 0.251n(^Ulnr 1TeV ^ 100GeV ^ /3 ' ^M a ,, 



ITeVV \ m d J\ T rh}2 J V ^1/12 



(5. 



While we expect this to be below the 60 e-foldings generally found for GUT-scale inflation, 
it need not be far smaller (e.g. by a factor of two), as one might have expected for a 
theory in which inflation happens at a much lower energy scale (e.g. ~ 1 TeV compared to 



9 We shall discuss the development of density perturbations during the "matter-dominated" phase 
of a Universe consisting of radions and other cold dark matter elsewhere, but note here that there 
is nothing special about the radion component, and it can be shown to behave as a typical dark 



matter particle. As shown in {TV, significant, additional fluctuations in the radion field are not 
produced during Phase II. 

10 This also guarantees that the radion energy density during cosmological nucleosynthesis was no 
more important than that of any other dark matter component, and therefore has negligible effect. 
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~ 10 12 TeV). This is because the Planck scale was relatively small during Phase I, when the 
main inflationary era occurred in our model. According to Eq. (|5.3p , it is not too difficult to 
satisfy the constraint that the number of e-foldings be considerably larger than one, unless 
o is absurdly small. 

Ours is only one of several proposed scenarios for cosmology in the brane world. The 
model expounded here has some overlap with that of Ref @, except that we assume that 
the radion and inflaton are different fields, resulting in substantial differences between the 
two models. By contrast to our model, Ref uses a bulk potential for the inflaton, leading 
to substantially different physics, and Ref ]T0| proposes a different theory of baryogenesis 
whereas we believe that baryogenesis from the (minimal supersymmetric) standard model 
electroweak phase transition is adequate. It should be clear to the readers that there may 
exist other viable cosmological scenarios in the brane world. Eventually, string theory should 
provide the appropriate radion/inflaton potentials, which hopefully will determine the cos- 
mological scenario that nature chooses. 

Recently, Randall and Sundrum proposed a scenario [22] where the extra dimension 
does not have to be compactified. In this scenario, the radius can have a run-away behavior. 
It is not even clear that the radion field has to be trapped during inflation to obtain the 
correct power spectrum of the density perturbation. The cosmology of such a scenario will 
be interesting to study. Some attempts along this direction can be found in Ref. p3| . 
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APPENDIX A: THE EFFECTIVE POTENTIAL 

Here we want to give some background on the motivation on the choice of the effective 
potential ( |1.2|) used in the text. Although our universe is non-supersymmetric, it is very 
helpful to start from a supersymmetric theory in which spontaneous supersymmetry breaking 
takes place dynamically. Since the brane world picture is naturally realized in string theory, 
where many non-trivial consistency properties (such as consistent quantum gravity) are 
automatically built in, we shall consider what stringy properties tell us about the effective 
potential. 

Since string theory has no free parameter (the string scale m s simply sets the mass 
scale), all physical parameters emerge as various scalar fields obtaining vacuum expectation 
values {vev ) determined by string dynamics. For example, the large radii of the large 
extra dimensions come from the vev s of the radion fields. Before supersymmetry breaking 
and dilaton stabilization (the latter fixes the string coupling value), the dilaton and the 
compactification radii are moduli, that is, the effective potential is flat (and remains zero) 
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as their vevs vary. This is true to all orders in the perturbation expansion. So supersymmetry 
breaking and moduli stabilizations are expected to come from non-perturbative dynamics, 
which is poorly understood at the moment. However, it is still reasonable to assume that 
the moduli degeneracy is lifted after dynamical supersymmetry breaking. 

Although the effective potential of a particular string vacuum (i.e., ground state) is 
model-dependent, there are stringy and supersymmetric features that are quite generic |24| . 
Here we shall give a very brief description of some of the properties that are relevant in 
this paper. Besides the graviton, the dilaton and the radii, a typical semi-realistic string 
model has gauge fields (in vector super- multiplets), charged matter fields (as components 
of chiral super-multiplets) as well as additional moduli. The general Lagrangian coupling 
N = 1 supergravity to gauge multiplets and chiral multiplets Zi (the index i labeling different 
chiral multiplets Zi will be suppressed) depends on three functions : 

(1) The Kahler potential K(z, z) which is a real function. It determines the kinetic terms of 
the chiral fields 



Ckin = K zz d^zd»z (Al) 

with K zz = d 2 K/dzdz. 

(2) The superpotential W(z) is a holomorphic function of the chiral multiplets (it does 
not depend on z). W(z) determines the Yukawa couplings as well as the F-term part of the 
scalar potential Vf : 

V F (z,z)=e K / M p^D z WK; 2 1 DW-3 l -^£^ , (A2) 

with D Z W ee dW/dz + WK z /M 2 Pl . 

(3) The gauge kinetic function f a b{z) is also holomorphic. It determines the gauge kinetic 
terms 

£„ = Ref ab F« u F^ b + lmf ab F^ b (A3) 
It also contributes to gaugino masses and the gauge part of the scalar potential Vd : 

V D = (Re/~ ^ (K z , T a z) (K- z , T b z) (A4) 
So the effective potential is given by 

V = V F + V D (A5) 

Consider a semi-realistic Type I string model, i.e., a D — 4, jV = 1 supersymmetric, chi- 
ral model, with a set of 9-branes and up to 3 sets of 5-branes, with a common 4-dimensional 
uncompactified spacetime (x to x 3 ). We shall treat the 6 compactified dimensions as com- 
posed of 3 (orbifolded) two-tori: the first torus with coordinates (x 8 ,x 9 ), the second with 
coordinates (xe,xr) and the third with coordinates (x^xs), the volumes of which are, crudely 
speaking, rf , r| and r| respectively. The 4-dimensional Planck mass Mp\ and the Newton's 
constant Gn are given by 

G N = M Pl (A6) 
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where A is the string coupling. The gauge couplings g$ and g$i of the gauge groups Gg and 
G 5i are 



_ 2 2{m? s rir2r^) 2 _ 2 2m 2 s rf 



g -2 = ^vp , ^ = ^p. (A7) 

where the ith set of 5-branes has as the size of its two compactified directions. For large 
radius r 3 , g 9 and g 53 are too small to be relevant, so the standard model gauge groups must 
come from the first two sets of 5-branes. This is the n =2 case. In the two examples PJ25 



that we know, this n =2 case is needed for phenomenology. Here, we shall identify r = r 3 
as the radion. To stabilize the moduli vevs (and maybe also to induce SUSY breaking), the 
string coupling A is likely to be strong. To obtain the weak standard model gauge couplings 
from a generic strong string coupling requires that m s ri and maybe m s r 2 to be around 10. 
This will modify Eq. (In semi- realistic string models, the picture is somewhat more 

complicated.) 

The dilaton and the volume moduli are bulk modes : 

S = gg 2 + id, T i = g^ + i6 i (A8) 

where the #'s are corresponding axionic fields. The radion field is parametrized by S and 
T 3 . For example, to lowest order, the gauge kinetic functions f 9 = S and f 5i = T i: while 
the Kahler potential is better known [24|. In the example of Ref H, with only two sets of 



5-branes (orthogonal to the third torus with very large r 3 = t~o), we have 



| z 12|2 



K = - ln(5 + 8- - £ n» ) - £ HTi + V) + ^-^L__ + ... (A 9) 

where z % i refers to open string chiral modes with one end of the string ending on the ith 
5-branes and the other end ending on the jth 5-branes. (For z u with i = 1,2, only the ith 
torus (world-sheet) excitation modes are included.) The superpotential W starts out with 
terms cubic in Zi 

W = y(S,T i ) jkl z :j z k z l + ... (A10) 

where y(S,Ti)jki are model-dependent functions of the moduli. 

Let us concentrate on the F term Vp of the effective potential. Generically, the lowest 
order terms in the brane mode effective potential are multiplied by some functions of the 
moduli, while higher order terms couple brane modes and the moduli. From the form of the 
superpotential W, where S and T 3 parameterize the radion field r, any brane potential will 
couple to the radion. In low orders, it will be a direct product of the brane potential and 
the radion potential. So the form V (ip)[l + fi(r)] is quite reasonable; this is the first term 
of our assumed effective potential (|1.2|). Choosing if) to be the electroweak Higgs field, the 
last term V\{ip) in Eq. ( |1.2|) is simply the Higgs potential in the standard model, at least 
for vev not much bigger than the electroweak scale. Of course, we are more interested in 
the electroweak Higgs potential in the minimal supersymmetric standard model, which have 
two Higgs doublets. There, the effective Higgs potential is only poorly known. 

Notice that Vp does not contain a term that involves only the moduli, which is a property 
that extends to all orders in the perturbative expansion. However, a term /o( r ) will appear if 
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some brane modes other than the inflaton develop non-zero vevs. Also, we do expect effective 
potential terms coupling the moduli to other bulk modes, as well as terms involving the 
moduli to be generated non-perturbatively. Otherwise, the moduli will appear as massless 
fields (much like the Brans-Dicke field), which is ruled out experimentally. Hence we need 
an effective potential to stabilize the radion. This is another reason we expect the presence 
of a term like fo(r). It is a bulk potential. This more or less justifies the choice of the form 
of the effective potential of the type ( |1.2| ) proposed in the text. 

APPENDIX B: BRANE INFLATION 

The brane inflationary scenario |jj emerges rather naturally in the generic brane world 
picture [0-0]. We may consider the Type I string where K branes sit more or less on top of 
an orientifold plane at the lowest energy state, resulting in zero cosmological constant. In 
cosmology, it is reasonable to assume that some of the branes were relatively displaced from 
the orientifold plane in the early universe. (This is the generic situation in F theory, which 
may be considered as a generalization of the Type I strings.) To simplify the problem, we 
assume that only one brane (or a set of branes) is displaced from the rest by a distance d. 
This situation probably arises after all except one brane have moved towards each other. 
Before supersymmetry breaking and dilaton stabilization, the force between the separated 
brane and the rest is precisely zero. In the realistic situation where supersymmetry is absent, 
we expect the potential V(d) to be, at large separation d, 

V(d) = m\d 2 ~ n {\ + e ~ <d ~ E e " mjd ) ( B1 ) 

where m[ are the masses of the NS-NS string states while rrij are the masses of the string 
RR fields (the sums are over infinite spectra). For large d and n — 2, V(d) is essentially a 
constant. The "1" term is due to gravitational interaction, the only long range force present 
at large d. For small d, the form of V(r) depends crucially on the mass spectrum. 

A key feature of brane inflation is the identification of the separation d with the vacuum 
expectation value of an appropriate Higgs field This Higgs field is an open string state 
with its two ends stuck on two separated branes. That is, this Higgs field is a brane mode 
playing the role of the inflaton. In the effective four- dimensional theory, the motion of the 
branes is described by this slowly-rolling scalar field, the inflaton -0 = m 2 d, which is the 
scalar component of one of the chiral field Zi, or some linear combination. To be specific, we 
shall at times consider the n = 2 case, and, as an illustration, keep only the graviton and 
one RR mode, resulting in an (over-) simplified effective potential, 

V(ip) ~ m*(l - e-W m ')F(if>) + \4(V> 2 ) (B2) 

where ipi and ip2 are two different brane modes and mj is a model-dependent mass scale, 
which is related to the mass of the RR mode m RR via ip\ = mid, that is, mjm RR = m\. We 
also include a generic smooth function F(i/j), which will be neglected in the text. Since the 
RR mode is massless before supersymmetry breaking, and that the supersymmetry breaking 
scale is below the string scale, we expect mi > m s > m RR . 

In this scenario, it is even possible that the electroweak Higgs field plays the role of the 
inflaton, a particularly interesting scenario. In this case, we may identify ip = ipi = ip2 as 
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the electroweak Higgs field and Vi(if)) as the electroweak Higgs effective potential. That is, 
Vi{ip = 0) ~ m% w . 



APPENDIX C: DERIVATION OF EFFECTIVE LOW-ENERGY DESCRIPTION 



In this appendix we derive the low energy, 4 dimensional description given in Sec. Li A 
above from the higher dimensional description of the brane world scenario. We start with 
the action 



L = J d s+1 x A ^-det(g AB ) 
where 



G 



K s 



+ J dV\/-det(^) \c h [ gflv (y x ), x(v X )} 



(s - 2)r( g /2) 

4( S - 1)VT S / 2 ' 



(CI) 



(C2) 



The notation here is as follows. The number of spatial dimensions is s = n + 3, n is the 
number of extra compactified dimensions, and G( s ) is the s-dimensional Newton's constant. 
The normalization of the first term in the action (|C1|) is chosen such that the force law 
at short distances is F = G( s )mim2/r s_1 . The quantities x A are coordinates in the higher 
dimensional space (the bulk) with < A < 3 + n, qab = 9ab{x c ) is the bulk metric, and 
is the Ricci scalar of #ab- In the second term, the quantities y M with < fi < 3 are 
coordinates on the brane. The induced metric on the brane is 



9^(V A ) =9A B [z c (y x )] 



dz A dz B 



(C3) 



where the location of the brane is x A = z A (y fl ). The quantity in (1) is the Lagrangian of 
all the fields, collectively called Xi th a t nve on the brane. 

This action ( |Cl"|) is a functional of the (4 + n) dimensional metric, of the location of the 
brane, and of the standard model fields ip, and is invariant under transformations of both 
the x A coordinates and the coordinates. We now specialize the x A coordinate system as 
follows. Lets write x A = (x^ 1 , x a ), where < /i < 3 and 4 < a < 3 + n. We can choose the 
coordinate system such that the brane location is 



^ = ^V) = (^,o,o,...,o). 



(C4) 



Hence we can identify the first four of the bulk coordinates x M with the brane coordinates 



We now make the ansatz for the metric, in the above choice of bulk coordinates x A , of 

(C5) 



g A B{x^,x a ) 



e 2 ^h ab (x c ) 



Here it is assumed that the internal space is a compact space of constant curvature, like an 
n-sphere S n or an n torus S 1 x . . . x S 1 , with metric h a iy(x c ). The volume V n and effective 
radius r of the extra dimensions are then given by 
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V n = r n = J d n x a sfdetihab). (C6) 
If r is the value of the radius r today, lets adopt the convention that 



d n x a ^/det(h ab ) = rl (C7) 
so that 

r = e*r . (C8) 

In going from the full metric gAB{x^i x a ) to the reduced form ( |U5| ) we have thrown away 
all the Kaluza-Klein modes which have masses > 1/r. Hence the ansatz ( |C5| ) will only be 
valid when all the fields vary with x M over length scales 3> r. We have also thrown away 
several of the components of the metric - the components g^ai^) and the traceless part of 
9ab{x^)- This is valid since these components have no couplings to the brane fields ip; in the 
four dimensional description they will act as free, massless scalar and vector fields which are 
coupled only to the metric g^ v Their equations of motion will be source free equations of 
the form V a V a ip — |26]], and so, at least classically, we can take them to vanish R. 



We now specialize the brane Lagrangian C b appearing in Eq. (|CT|) to be of the form 

X) = + £ rest (iV, Xrest), (C9) 

where ip is the inflaton field or fields, and x rest denotes the remaining brane fields other than 
if), described by the Lagrangian C Tes t- We also add to the action the terms 

- J d s+1 x A ^-det(g AB ) Kuik($) - / dV^-det^) V hlwe (iJ>, $). (CIO) 

This consists of a bulk potential energy per unit s-dimensional volume Vbuik an d a brane 
potential energy per unit 3- volume Vb ran e Q- ^ e discussed in Appendix |A] above the physical 



11 To see that these fields are exactly decoupled classically, one should use a definition of the radion 
field $ which is more general than Eq. (|C5|), namely 

J d n x a ^/-det(g AB ) = r"e n *^-det( 5 ^). 

This equation together with the action (CI) shows that $ is the only piece of the metric which has 
couplings to anything other than the metric. 

12 Quantum mechanically, these fields will be subject to the same process of parametric amplifi- 
cation during inflation as normal gravitons, and if they start in their vacuum states the total in 
these fields today should presumably be comparable to the total O in relic gravitons from inflation, 
which is of the order of 10 -14 in typical inflation models but smaller in the models of this paper. 

13 Note that the explicit dependence of these potentials on the metric component $ spoils the 
covariance of the full action under transformation of the x A coordinates; it is difficult to write 
down a fully covariant radius-stabilization mechanism. 
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origin for such terms which depend on the size of the extra dimensions as well as on the 
infiaton. 

Using the ansatz (|C5|) in the action (|Cl|), inserting the brane action ( |C9| ) and adding the 
potential terms ( |C10|) now yields the reduced 4-dimensional action 

S = J d'x^Tg {^ {A) R + "^^^^y - V{1>, *) + CUSw 4>)] ■ (Gil) 



Here G is the usual 3-dimensional Newton's constant, given by 

1 K„rJ} 



(C12) 



167TG G {s) ' 

where tq is the equilibrium value of the radius of the extra dimensions. The action ( |G11|) 
has the form of a scalar-tensor theory of gravity, written in the Jordan frame. Note that the 
sign of the kinetic term for the scalar field in the action ( pll|) is opposite to the normal sign; 
this is not a problem since it is the sign of the kinetic energy term in the Einstein frame 
(see below) that is relevant to considerations like stability and positivity of energy etc. The 
Jordan-frame potential V is given by 

V(il>, <J) = r V^ bu ik($) " 7^r- 2 e (n " 2)l> + V biane (^, $), (C13) 

07TLtTq 

where the Ricci scalar of the metric h a b is and ki is a dimensionless constant of order 

unity (cf Eq. ( |C7| ) above). From now on we specialize to flat internal spaces so that ki = 0. 
Then we see that only the particular combination V = r n Vb u \k + Hrane of the potentials Vbuik 
and Vbrane is relevant in the low energy description. Our assumed form for this potential V 
is given in Eq. ( |1.2|) above. 

Finally, we transform to the Einstein frame description. We introduce a canonically 
normalized radion field by defining fi according to Eq. (|2.3|) above, and we define 



$ = n^$. (C14) 

The radius r is thus related to <l> by 



r exp 



(C15) 



from Eq. ( |C8|) . We define the Einstein frame metric by 

gap = ^^g aP = e ni g aP . (C16) 
The action then takes the form of Eq. ( |2.1| ) above, where 

<S , rest[fl'a/3,Xrest] = / ^x/-^Aest(^o/3> Xrest) (C17) 

is the action of the brane fields Xrcst- 

The equations of motion derived from the action Q2.1|) , when we treat the last term as a 
fluid with Jordan-frame density p and pressure p, are 
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M 2 ~ 1 

S7T Z 



— e 



-2$ 



+ e" 2 ^ [(p + p)^ + pg af} ] , 



(C18) 



V a V Q $ + i-e-^(W) 2 - A [e" 2 */^(^, $)] + -U^Go - 3p) = 0, (C19) 



and 



V a V> - -V a $ V> - e-^^-rVN), $) = 0. 
p o*0 



(C20) 



Here V a is the derivative operator associated with the Einstein frame metric ^ a/ g, and u a is 
normalized with respect to g a p. 



APPENDIX D: A PARTICULAR REALIZATION OF PHASE EVOLUTION 

Here, we study the pre-inflation Phase in some detail for the sinusoidal potential 

fi(r) = a(l — cos mr) (Dl) 



introduced in § [II B| . For this choice, the effective potential for the radion field $ in the 
Einstein frame is, from Eq. (|2.14j) , 



r 



2n 



Kfr = Voy — j (1 + a — acosmr). 

If amr/2n ^> 1, it is easy to show that the maxima of V e s are at 

2n{l + 2a) 



(D2) 



mrj = (2j + 1)tt + 6j 



amrj 



(D3) 



The heights of successive maxima of V e g differ by approximately 

'r Q \ 2n (r j+l - rA 



Vtl +1) -V$ ^-2nV {l + 2a) 



(D4) 



the potential is biased toward large values of r. Since whether or not the radion escapes over 
the first barrier is determined at very early times, we can ignore expansion, and treat the 
dynamics as completely conservative. If the radion starts out with zero or negligible kinetic 
energy, then if it can only escape over the barrier at r J+1 if it begins sufficiently close to rj. 
Expanding the effective potential near rj we find that 



V eS (rj + Arj) - y$ « -2naV ( — ) (mAr J 



r 



2/i 



(D5) 



V eS (rj + Arj) = V eff (r i+ i) when 
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mArj ~ 



\ 



8t™(1 + 2a) 

(D6) 



amrj 



which is small if amrj/2n ^> 1. If the radion begins at any radii rj + Ar,- < r < r J+ i it 
should be captured at the nearest minimum, which is roughly halfway between rj and r^+i; 
if it begins at rj < r < rj + Ar^, then it should grow without bound. Numerical solutions of 
eqs. ( ^.17| ) verify this picture. If the value of r at the beginning of Phase is random, the 
probability that the radion is not trapped is P e ~ mArj/2n, which decreases oc (2n/ amr) 1 ^ 2 
as amr/2n grows. Thus, for large amr/2n, it is extremely likely, although not guaranteed, 
that the radion is trapped at the minimum of V e g nearest to its starting value of r. 

The Einstein- frame expansion rate once the radion settles into a minimum of V e s(Q) is 

^ 2 -^exp(-2$/^), (D7) 

where $ is the value of the field. Expanding V c fi(<&) around the minimum implies 

- (mr\ 2 (A$) 2 
V ee $ + A$ )»aVb — exp(-2$/(i)±—l-, (D8) 
\riji J 2 

so the characteristic oscillation frequency for fluctuations in the radion field is 

9 fmr\ 2 , . i 12a(mr) 2 H 2 , 
wjwo — exp -2$//i « 1 / . (D9) 

For large values of mr (and a not too small), c<j$ ^> H . One consequence of this inequality 
is that we do not expect large-scale fluctuations in the radion field to arise during inflation. 
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